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ABSTRACT 

Weyl focusing effects on the image magnifications are investigated by using the mul- 
tiple gravitational lens theory. We focus on the gravitational lensing effects due to the 
small scale virialized objects, such as galaxies and clusters of galaxies. We consider a 
simple model of an inhomogcncous universe. The matter distribution in the universe is 
modeled by randomly distributed isothermal objects. We found that, for the majority 
of the random lines of sight, the Weyl focusing has no significant effect and the image 
magnification of a point like source within redshift of 5 is dominated by the Ricci 
focusing. 
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1 INTRODUCTION 

Distance-redshift relation plays an important role in the astronomy and observational cosmology. The standard distance, 
which has been used in the most of previous studies, is based on a postulate that a distribution of matter in the universe 
is homogeneous (e.g., Weinberg 1972). It has been, however, well recognized that our universe is highly inhomogeneous on 
small scales. Since the inhomogeneities of the mass distribution focus (defocus) the bundle of light rays (the gravitational 
lensing effects), the distance in an inhomogeneous universe deviates from that in the homogeneous Friedmann universe. It is, 
therefore, obvious that a detail understanding of the propagation of light rays in the inhomogeneous universe is necessary for 
correct studies of objects in a distant universe. 

Since the pioneering work by Gunn (1967), there has been a lot of progress in this subject. Babul & Lee (1991), among 
others, studied the effects of the Ricci focusing by weak inhomogeneities. They found that the dispersion in image magnifica- 
tions due to large scale (<; 0.5/i _1 Mpc, where Hq = 100/ikm/sec/Mpc) structures is negligible even for sources at redshift of 4. 
They also pointed out that the dispersion is very sensitive to the nature of the matter distribution on small scales. The same 
result was also obtained by Frieman (1996). He improved Babul & Lee's study to reflect the recent developments in numerical 
and observational studies of the large scale structure. In these studies, the effects of the Weyl focusing which induce a shear 
of light ray bundle were neglected. Nakamura (1997) examined the effects of the shear on the image magnification in the cold 
dark matter model universe with linear density perturbation. He found that the effect is sufficiently small and concluded that 
the Weyl focusing can be safely neglected for a light ray passing through a linear density inhomogeneities. Jaroszyriski et al. 
(1990) and Wambsganss et al. (1997) used the multiple gravitational lens theory with a large N-body simulation of the cold 
dark matter universe. Their result agreed well with analytical studies. 

The above studies mainly focused on the large scale inhomogeneities, whereas the effects of small scale objects, such as 
galaxies and clusters of galaxies, have not been fairly taken into account. Kayser & Refsdal (1988) investigated the gravitational 
lensing effects due to randomly distributed King model galaxies. They paid a special attention to a high magnification part 
of the magnification probability distribution. Recently Wambsganss, Cen & Ostriker (1998) studied the gravitational lensing 
effects by using the large N-body simulation with an effective resolution of comoving 10ft~ 1 kpc. They, first, shoot the light ray 
through the lens planes by using the multiple gravitational lens equation, then the magnification matrix is determined from 
the mapping of the light ray positions between the image and source plane. In this procedure, the Ricci and Weyl focusing 
can not be treated independently, therefore no discussion is given for the Weyl focusing effect in their paper. However, the 
magnification factor as a function of position in the source plane and image plane are presented in Figure 4 and 6 of their 
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paper. Those Figures show that there is quite a large region that is demagnified by a small amount, and a few relatively small 
spots that are quite highly magnified by the small scale inhomogeneities. 

The purpose of this paper is to examine the Weyl focusing effect due to the small scale inhomogeneities, and we are not 
concerned with the effect due to the large scale structure. The density distribution in the universe is modeled as a randomly 
distributed isothermal lenses. The isothermal lens is approximated by virialized objects, such as galaxy and cluster of galaxies. 
This model is similar to the one studied by Kayser & Refsdal (1988). However our emphasis is different from theirs, i.e. we 
pay special attention to the gravitational lensing effects on a majority of light rays. Although this model is a very simplified 
and unrealistic one, we believe that the model is good enough to investigate the essential points of the Weyl focusing effect 
due to the small scale inhomogeneities. 

This paper is organized as follows: fn section 2, we describe our models and method of simulating the light propagation. 
The results of the simulation are summarized in section 3. The paper concludes with discussion in section 4. 



2 MODELS AND METHOD 

2.1 Theory of multiple gravitational lensing 

We use the multiple lens equations to trace the propagation of infinitesimal bundles of light rays. Schneider, Ehlers & Falco 
(1992) deals with the theory of multiple gravitational lensing in detail. Here we simply describe only the aspects which are 
directly relevant to this paper. 

As was done in previous studies (see e.g., Blandford & Narayan, 1986 and Kovner 1987), we consider N screens (lens 
planes) between an observer (z = 0) and source (z = z 3 ), located at redshifts z< with i runs from 1 to N + 1 with zjv+i = e s . 
In the following, the quantities on the lens and source planes are described by indices {A, B, ..} = {1, 2}. The position vector 
of a light ray on the i-th lens plane is denoted by yA(zi). Let &A(y(zi)) denotes the deflection angle of a light ray at position 
yA(zi) on the i-th lens plane. The multiple gravitational lens equation and the evolution equation of the lensing magnification 
matrix Mab are written as, 

2- 1 

VA(zj) = -j±yA(zx) - ^ Dija A (y(zi)), (1) 

3— 1 

MAB{y{ Zj )) ^Sab-Y, ^^&A,c(v{z i ))McB(y(z i )), (2) 

i — l ^ 

for 2 < j < N + 1, where Dij (Di) denotes the standard angular diameter distance between redshifts of Zi and Zj (0 and Zi) 
with i < j, the comma denotes differentiation with respect to the components of yA(zi) and y(zi) = yA(zi). The deflection 
angle &A(y) is determined by the following equation: 



& A (y) = %{ dV J Mt2 W) " <E» » (3) 
J \y-y\ 

where E(y) is the surface mass density and (E) is its average value. For a notation convenience, we introduce the following 
quantities: 

a A {zi,Zj) = Dij&Aiyizi)), (4) 
T A B(zi,Zj) = D p Dl & A ,B(y(zi)). (5) 

The optical tidal matrix Tab(zi, Zj) is decomposed into the Ricci and Weyl focusing terms, respectively: 

K{zi,Zi) = ^{T 1 (z 1 ,z,)+T 22 (z t ,z J )), (6) 

T{zi,zj) = - (Tn(zi,Zj) -?ii{zi,Zj)) + Wi2{zi,Zj). (7) 

In general, equation (^) is not an explicit equation for Mab, since the equation involves a summation over Tab evaluated 
on the light ray path, such that one first has to solve the multiple gravitational lens equation (|l|). However for the light rays 
traveling in regions where a a and Tab < 1, one can expand Mab in powers of a a and Tab about its value when the light 
ray is unperturbed. We rewrite equation ([j]) as yA(,Zj) = Ha' + Va^( z 3) + ^i - 2 ) an d equation (^) as M.AB{y{z 3 )) = 
M { AB(y (0) fe)) + M%(y m (Zj)) + MA%{y m (zj)) + where y { 2\zj) is the first term of the right hand side of equation 
(0) and y^P(zj) is the second term, but the deflection angle is evaluated at the unperturbed light ray position. Expanding 
equation (|2|) in terms of a a and Tab, one finds 

Mf B { Zj )=5 A B, (8) 
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3-1 i-1 



Mab( z j) - E 5^ [7Ao(«i, Zj)TcB(z k , Zi) + T A B,c{zt, Zj)ac(z k , Zi)) , 



(9) 



(10) 



for 3 < j < iV + 1. In the above expressions, 7as and a a are evaluated at the unperturbed light ray position. The image 
magnification factor of a point like source is given by the inverse of the determinant of the magnification matrix, i.e. p, = 
|det Mab I ~ X • Up to the order of M-ab, the determinant is 



det TWas(zj) = 1 - 2^2'R(zi,z j ) + 



i=l 
3-1 i-1 



3-1 



E 7 ^ 



3-1 



+2 2_, /J {TZ(zi, Zj)TZ(z k ,Zi) + Re[.F*(zj, z j )J r (z k , Zi) ]+Tl,A(zi, Zj)aA(z k , Zi)} . 



(11) 



for 3 < j < N + 1. This is our principal equation. In general, the effects of Ricci and Weyl focusing on image magnification are 
coupled. However up to this order, they are not coupled. We call the terms in the equation ( |ll| ) which involve Ricci focusing 
terms as "Ricci contribution" and that involve Weyl focusing terms as "Weyl contribution" . 



2.2 Truncated singular isothermal sphere lens model 

We adopt the truncated singular isothermal sphere as the lens model. Its surface mass density is written as 



^ = ( 



R Y 



(12) 



where a v is the one-dimensional velocity dispersion and Rg is the half mass radius (Pei, 1993). The corresponding 3-dimensional 
mass density runs as r~ 2 (r Rg) and as r~ 4 (r 3> Rg), and the mass within the radius R is 

. , TvalRc R 
m ^ R)= G Rc + R > 

then the total mass is m to t = ttu^Rg/G. This models should provide a fair approximation to virialized objects, such as a 
galaxy and cluster of galaxies with isothermal dark halos. 

For the distribution of matter in the universe, we assume that the isothermal lenses are randomly distributed with the 
average mass density phiz) and the rest of the matter has the uniform distribution. We also assume that the comoving density 
of lenses is constant in time, thus ph{z) = (1 + z) 3 pl(0). Furthermore, we approximate the lensing effects of the isothermal 
objects, except the nearest one, to a form of uniform surface mass density. 

Under the above assumptions and the circularly symmetric mass distribution of the lens model, the deflection angle is 
given by 



a{zi,Zj) = \Ja\{zi,Zj) +al(zi,Zj) 



An ( — 



Di 



Rg 



R 2 c (Rg+Rc 



Rg + R 

Rg 

Rg + R R c (Rg + Rc 



R G R 

Rc 



R g R 



with 

O'er (Zi , Zj ) 



47T 



(t) : 



(13) 



(14) 



where 7? is the distance between a light ray position and center of the nearest lens object in the i-th lens plane. R c is a half 
of the mean separation length between the lens objects defined by 



Rc 



mtot 



ttY,] 



= \J a cr (zi, zj)Rg 



-n L (o)(i 



2 DjjDj cdt 
Di dz 



(15) 



where is the average surface mass density of the lens objects in the i-th lens plane, and 5 z is the redshift interval between 
the (i — l)-th and i-th lens planes, and Ql{0) = pi(0)/p cr (0) = pl(0)(8tvG /SHq) is the density parameter of lens objects. 
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The first term in square bracket of the second line in equation ( |13[ ) describes the deflection due to the nearest lens object, the 
second term is due to the others, and the last term is due to the background average mass density of lens objects. From the 
deflection angle (|l^), the Ricci and Weyl focusing terms are immediately given by 



K(zi,Zj) = 



■■{Zi 



R 



Rg 



2 R{R G + R) 2 
1Rg(Rg + 2R) 



Rc(Rg + Rc 



2 R(R G +R) 2 
We introduce a compactness parameter v as follows; 



(16) 
(17) 



{Oy\ 2 1 C 



( lOkpc 



(18) 



^200km/sec / 

This parameter measures the effectiveness of an isothermal object as a gravitational lens. As was pointed out by Kayser & 
Refsdal (1988), the model with randomly distributed isothermal lens objects is completely described by two parameters. The 
parameters they used depend on the distance between source and observer. On the other hand, we characterize the model by 
the compactness parameter v and density parameter of lens object f2i(0) (hereafter we shall simply denote it as Ql) which 
are independent of redshifts of the source and lenses. 

We should note that the optical tidal matrix Tab becomes larger than unity at the very central region of the lens object 
(the inner region of ~ 0.5 x Einstein radius), therefore the validity of assumptions used in deriving the equation (jnj) breaks 
down for light rays passing through that region. We now estimate the strong lensing effect on our study in two ways. 

The first is based on an order-of-magnitude estimate (section 14 of Peebles, 1993, and also Futamase & Sasaki, 1989), 
we examine a magnitude of the tidal matrix Tab- Suppose the lensing objects are randomly distributed and each with 
mass M = 2a v 2 l/G, where I is a characteristic comoving size of a lens object and is of order Rg- Hence the mean co- 
moving number density of the lens objects is til = fiL(3-ffo 2 /167r)cr„ _2 Z _1 , so the mean comoving separation distance 
is ro = S!L~ 1 / 3 (3i/o 2 /167r)~ 1/ ' 3 <T t , 2 ' /3 Z 1 ' /3 . Then for a geodesic afhne comoving distance of A, the light ray gravitation- 
ally encounters such objects N g = A/ro times in average. At each encounter, the contribution to the tidal matrix is 
5T = 4Tv{a v 2 /c 2 ){r /b 2 ){D d D ds /D s ) ~ 4n{a v 2 /c 2 )^ 1 {D d D ds / D s ), where D ij is a comoving angular diameter distance 
with the subscript d (s) stands for a lens (source), b is the comoving impact parameter, and we have assumed that the mean 
comoving impact parameter is of order ro. Since the sign of each contribution will be random, the total contribution to the 
optical tidal matrix will be 



3/4VTTc 



4tt(- 



,al c 



i Ho D d D d 



H 



Ho D d D di 



^A 



(19) 



The contribution from the direct encounters can be similarly estimated by noting that the average number of encounters is 
Nd = / 2 A/ro 3 = \Ql{3Ho 2 /16tv)o- v ~ 2 1, with each encounter contributing 5T d — 4n(a v 2 /c 2 )l~ 1 (D d D ds /D s ) with random sign. 
The result turns out to be the same as that of gravitational distant encounters, equation (|l9|). In the case of Einstein-de 
Sitter universe model, the comoving distance A becomes c/Ho at the source redshift z s — 3 and the averaged value of distance 
combination over the lens redshifts is {(Ho/c)(D d D ds /D a )) = 1/6 for any redshifts of source. Accordingly, we find that the 
magnitude of the optical tidal matrix scales as ~ 0.2V^l\/^ for the source redshift of 3. In the following we only consider 
lens models with Q.l < 1 and v < 1, therefore a typical value of the optical tidal matrix can be expected to be of order O(0. 1) 
or lower. Thus, we can conclude that Tab is less than unity for a majority of random lines of sight. 

Next we examine the lensing optical depth defined by Turner, Ostriker & Gott, (1984) to quantify the probability of the 
light rays being affected by strong lensing. The differential optical depth for the truncated singular isothermal sphere model 
is given by 



dr 



-n L E(i + z d ) 



D d D ds cdt 
D s dz d 



dz d 



where 



E = 



IRg 

4 o cr 



+ 



- 1 



(.'- 



R C (R G +Rc) 



(20) 



(21) 



We numerically integrate the last equation for cases of the lens models with (Ql,v) — (1,0.1), (1, 1), (0.2,0.1) and (0.2, 1), 
and for the Einstein-de Sitter universe model. The results are presented in Table 1. From the Table 1, it can be found that the 
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Table 1. The lensing optical depth. 



2 
3 
1 

5 



= 1 



(«£,!/) = (1,0.1) 



9.8 x 10~ 4 

2.9 x 1CT 3 
4.9 x 1CT 3 
6.6 x 10" 3 

8.1 x icr 3 



(1.1) 

8.5 x icr 3 
2.5 x 10" 2 
4.1 x 10~ 2 
5.5 x icr 2 
6.7 x 10~ 2 



(0.2,0.1) 
2.0 x 10~ 4 
5.9 x 10~ 4 
9.8 x 10" 4 
1.3 x 10" 3 
1.6 x 10" 3 



(0.2,1) 
1.7 x 10~ 3 

5.0 x 10~ 3 

8.1 x 10" 3 
1.1 x 10~ 2 
1.3 x 10~ 2 



probabilities of strong lensing events are very small except for an extreme model with (Ql, v) = (1, 1). Even for the extreme 
model, the probability is not significantly large. We thus conclude that the strong lensing effects do not significantly alter our 
results, ft can be said from the above two estimations that we can safely use the perturbative equation (JTTJ) . 

2.3 Ray shooting 

Since we have assumed the random distribution for lens objects, the probability of finding lenses in some region on a lens plane 
is described by Poisson distribution. If one sets the redshift interval 5z to be sufficiently small, the surface number density 
of lenses becomes small. In this case, the lensing effects are mainly due to the nearest lens and are well approximated by the 
equations (|l3]), (|lrj ) and (|f7|). Then the all necessary information about evaluating the magnification factor ( |ll| ) is obtained 
by randomly determining the relative position between a light ray and the nearest lens object in each lens plane. We perform 
Monte-Carlo simulations to trace the propagation of light rays. The procedure is described in the following: 

(i) First of all, we determine the redshift intervals of lens planes to satisfy the condition that Rg /Rc is sufficiently small. 
We set Ra/Rc < 0.1. 

(ii) The relative positions between a light ray and a center of lens object are randomly determined in each lens plane. 

(iii) Summations in equation ([y]) are performed for each term, and the results are stored in a file. 

Steps (ii) and (iii) are repeated for each light ray. 
3 RESULTS 

For the background universe, we only consider the Einstein-de Sitter universe model, i.e., Q.o = 1 and Ao = 0. For the lens 
models, we choose the compactness parameter v = 1 and 0.1 which roughly correspond to a galaxy scale inhomogeneity and 
the scale of a cluster of galaxies respectively. The density parameter of the lens objects are set to be Ql — 1 and 0.2. We 
choose the source redshifts z s — 1, 2, 3, 4 and 5. From the condition Rq/Rc < 0.1, the redshift intervals between lens planes 
are typically set to be ~ 10~ 2 . 

10 6 runs are performed for each model. For each run (light ray), we then obtain the Ricci and Weyl contribution and 
identically the image magnification factor. This immediately gives distribution functions of runs on the Ricci- Weyl contribution 
plane for all the source redshifts. We calculate number densities of results of runs in Ricci- Weyl contribution plane. The peaks 
of the number density and the isodensity contours which enclose 68% of all runs are presented in Figure 1. The probability 
distributions of Ricci and Weyl contribution for the case of the source redshift of z s = 3 in (Ql, v) = (1, 1) model are shown 
in Figure 2. As is clearly shown in the optical scalar equation (see e.g., Schneider et al. 1992), the Weyl contribution is always 
negative. Figure 1 reveals that the Weyl contributions in the majority of light rays are rather small except in the case where 
the source redshifts z 3 > 3 in v) = (1, 1) model. Alternatively, this point can also be shown in Figure 2. The probability 
distribution of the Weyl contribution has a narrow peak centered at very small values, in marked contrast with that of the 
Ricci contribution which has a broad distribution. Since we have assumed no evolution for lens objects, the dispersion keeps 
on spreading in Ricci- Weyl contribution plane even for high redshift. 

In order to examine the effects of the Weyl focusing on the image magnification quantitatively, we calculate the magni- 
fication factors evaluated without the Weyl contribution (denoted by /j.r). Then we calculate the following quantity which 
measures the influence of Weyl contribution on the image magnifications: 



where the summation is taken only over rays with det Mab > 0. The rays with det Mab < belong to a multiple image 
system, therefore the percentage of runs which are excluded in the above evaluation roughly exhibits the probability of the 
strong lensing events among random lines of sight. At the same time, the above condition also excludes the rays which pass 
through the high Tab region. The results are presented in Table 2 with the percentages of the excluded runs (in parentheses). 
Observationally, the strong lensing events among the random lines of sight are very rare. For example, the probability of 




(22) 
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Figure 1. Isodensity contours which enclose 68% of all runs and peaks of surface number density of runs. The peaks are denoted by 
pluses. The solid lines are for a source redshift z s = 1, dashed lines for z s = 3 and long-dashed lines for z s = 5. The lens model parameters 
(O^; u) are denoted in each figure. 



Table 2. Afi, the effect of the Weyl contribution on the image magnification with the percentage of runs excluded in this evaluation (in 
parentheses) . 

{n L ,v) = (1,0.1) (1,1) (0.2,0.1) (0.2,1) 

z 3 = 1 1.8 x 10~ 2 (0.16) 5.7 x 10" 2 (1.0) 7.6 x 10" 3 (0.024) 2.6 x 10" 2 (0.18) 

2 3.1 x 10- 2 (0.55) 1.0 x 10" 1 (3.3) 1.4 x 10" 2 (0.083) 4.3 x 10~ 2 (0.56) 

3 4.0 x 10~ 2 (0.96) 1.3 x lO^ 1 (5.6) 1.8 x 10" 2 (0.15) 5.5 x 10~ 2 (0.95) 

4 4.8 x 10~ 2 (1.34) 1.5 x 10" 1 (7.5) 2.1 x 10~ 2 (0.21) 6.5 x 10" 2 (1.3) 

5 5.5 x 10~ 2 (1.7) 1.7 x 10" 1 (9.2) 2.3 x 10" 2 (0.27) 7.2 x 10" 2 (1.6) 

the multiply imaged quasars in a quasar sample is at most 10 -2 (e.g., Claeskens, Jaunsen & Surdej, 1996). Combining this 
fact with our results summarized in Table 2, it may be reasonably concluded that, as far as our simple matter distribution 
model is concerned, the Weyl focusing has a negligible effect on the majority of light rays even for sources at the redshift of 
5. This result can be naturally explained by the following two reasons: First, the Weyl focusing is a second order effect on the 
image magnification, therefore it becomes important only for the light rays passing through a very high non-linear (relatively 
rare) region. Secondly, since we have assumed the random distribution for the isothermal lenses, the rays coherently affected 
by the Weyl focusing are very rare, consequently, the majority of light rays are only weakly influenced by the Weyl focusing. 



4 DISCUSSIONS 

In this paper, we restricted our study to the gravitational lensing effects due to the randomly distributed isothermal lenses. 
We found that the Weyl focusing effect is small, A^i ^ 0.1, for a majority of light rays. 

Lee & Paczyhski (1990) examined the gravitational lensing effects in randomly distributed clumps with Gaussian surface 
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Figure 2. Probability distributions for a model with 
= (1; 1) and z s = 3. (a) Isodensity contours 
which enclose 38%, 68% and 87% (inner to outer) of 
all runs. The plus denotes the peak of the surface num- 
ber density of runs, (b) The probability distribution of 
Ricci contributions, (c) The probability distribution of 
Weyl contributions. 



mass density profile. They found that the image magnifications are dominated by Ricci focusing and the Weyl focusing has no 
significant effect. They only considered a case of a source redshift of z s — 1.631, in an fi u = 1 and An = universe. However 
we have found that, as far as the random distribution of lens objects is concerned, the Weyl focusing is rather small even for 
higher redshift. 

In this study, a correlation of lens objects and large scale structure are not taken into account. The study of the influences 
of the correlation on the Weyl focusing lies outside the scope of this paper, and will be examined in future works. On the other 
hand, the Weyl focusing effect due to the large scale structures is investigated by using N-body simulation (Jaroszyriski et al., 
1990) and analytically (Nakamura, 1997). These two studies show that, although there is an uncertainty of a normalization 
of the density power spectrum, the magnitude of Weyl focusing due to the large scale structure is of the order 10 -2 ~ 1CP 1 
for the source redshifts of 1 < z s < 5 (Figure 3 of Jaroszyhski et al., 1990 and Figure 3 of Nakamura, 1997). Consequently it 
can be said that, comparing the results of the above mentioned studies with our results, the Weyl focusing effect due to the 
large scale structures is comparable with or larger than that due to the small scale inhomogeneities. 
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